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La métrologie existe depuis plus d'un siècle mais, à quelques 
exceptions près, cette discipline ne s'est appliquée à la chimie 
que durant les quinze dernières années. 

Ce module permet à l’étudiant : 

4 - de savoir identifier et affermir la place de la fonction 

"Métrologie" au sein d’un laboratoire de chimie, en relation 
avec le système d'assurance qualité en vigueur ou en projet 
dans cette structure, 

4 d’apprendre à avoir confiance et inspirer confiance dans 
des résultats de mesure, d'analyse ou d'essais, 

’i-de maîtriser les outils associés, 

'4-d’acquérir les fondements de base sur le contrôle de qualité 
et sa mise en application, 

4 - de compléter sa formation dans les techniques 
d’expression et de communication- 
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El : Métrologie E2 : Assurance qualité 


I. Définitions 
IL Précision et exactitude 

III. Introduction à la pratique des 
Plans d’expériences 

IV. Plans de pesée 

V. Matrice d’Hadamard 

VI. Plans factoriels complets 


I. Démarche qualité 

II. Carte de contrôle 

III. Justesse 

IV. Répétabilité 

V. Test de linéarité 
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El 



Historique de la métrologie 


"L'homme est la mesure de toutes choses" 
Protagoras (sophiste grec, 485-411 av JC) 


"La métrologie est une des bases de la qualité" 




^Les unités de référence : 

Besoin d’unités de référence (unités de temps, poids, volume, température...). 

Une classification arbitraire des quantités fondamentales de longueur, masse et 
temps qui ont déjà une longue histoire. 

- Exemple : l’origine de la mesure du Yard anglais était la longueur du bras du roi de 
l’époque (Henri 1 er ; 1100-1137). Il correspondait à la distance entre le bout des doigts 
et le nez du roi. 

-Le premier « Yard étalon », constitué d’une barre de fer, a été établi entre (1272- 
1307). D’autres Yards étalons constitués (en argent ou en bronze) au cours des siècles 
suivants. 

-Le Yard Impérial a été établi en 1878 par l’acte sur les Poids et Mesures avec des 
spécifications précises. 

Des bases de mesure : 

- Le pied de roi (o, 32483 m ) se subdivise en 12 pouces, 

- le pouce ( 2,706 cm) se subdivise en 12 lignes, 

- la ligne ( o,226 cm ) se subdivise en 1 2 points 

- le point ( o,i88 mm ). 
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• Quelques membres de l'Académie Royale des 
Sciences ralliés à la révolution, proposent 
l'établissement d'une unité de longueur 
universelle sur des base géodésiques vue les 
différences de mesures d'une région à l’autre. 


•Le 26 mars 1791 naît le mètre , dont la longueur 
est établie comme égale à la dix millionième 
partie du quart du méridien terrestre .(la 
circonférence de la terre est rigoureusement 
égale à 40 000 km). 
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Métrologie Chimique: Validation des 
Méthodes d'Analyses Chimiques 

Métrologie en chimie est devenu une référence utilisée par les 
techniciens et chercheurs de laboratoires d'analyses, et d'instituts de 
recherche. Il est aussi particulièrement recommandé aux ingénieurs et 
étudiants de 3^ cycle (mastère, doctorat) en chimie analytique, 
environnement et qualité. 

•Objectifs: 

Acquérir les outils nécessaires pour caractériser une méthode (tests 
statistiques, planification des expériences...). 

Savoir organiser les essais pour les optimiser d'un point de vue 
technique, mais également économique. 

Utiliser les données issues des travaux de caractérisation pour 
estimer l'incertitude du résultat de l'analyse et valider cette méthode. 

•Prérequis: 

Connaissances de base en statistiques souhaitées. 
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Chapitre I 
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1/ Définition de la Métrologie 

>^La métrologie peut se définir comme étant "la science de la mesure associée à 
l’évaluation de son incertitude”. Au sens large; c’est une science de la mesure . 



Mesure = résultat d 

a 

u mesurage d’un mesurande 


Résultat de la mesure Action d< 

mesurer Grandeur n 

Mesurée 


S Ce terme est utilisé dans le domaine des mesures de paramètres physiques (par 
exemple température, pression, dimensions, temps, etc.), et plus fréquemment 
employé dans des calculs de constantes d’équilibre chimiques 
(identification/quantification d’éléments ou composés) ou biologiques (par exemple 
comptage d’organismes microbiens). 

S C’est une composante de la qualité. 

^Sa spécificité est dans la validation du résultat. 

^La chimie analytique représente une des pierres de base de la chimie et englobe le 
concept récent de « métrologie chimique » dont les principes ont été trop souvent 
sous-estimés par les analystes. 


Objectif : Obtenir des mesures f^l^^-^eau de confiance élevé) 
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La métrologie présente trois aspects 


/ 

\ 

Métrologie des équipements 

Contrôle qualité des mesures 

> Vérification du fonctionnement 
correct des appareils par mesure de 
grandeurs fondamentales 

(ex : température, pression, masse, volume. . .) 

^Nécessité 

■locaux adaptés, 

■personnel formé et compétent 

> Détermination 
des incertitudes 

sur les mesures -> sur 
les analyses 

Validation des méthodes 
« alternatives » 


La métrologie et ses exigences sont intégrées dans les normes qualité 
concernant les laboratoires : 1 7025 et BPL. n 



1.1/ Métrologie des équipements 


^Vérification du bon fonctionnement des équipements avec : 

- Rendus de mesures fidèles et justes 

- Respect des spécifications du constructeur. . . 

-Mesures fidèles : même résultat obtenu ou mesures voisines en répétant le 

mesurage plusieurs fois. 

-Mesures justes : mesure en accord avec la valeur attendue. 

^Vérification de chaque type d’appareil ( balance, thermomètre, 
spectrophotomètre...) selon un cahier des charges : 

- extrêmement précis et rigoureux, 

- faisant l’objet d’une norme constructeur. 

Un cahier de vie accompagne chaque instrument 


Document qui recense toutes caractéristiques et les interventions faites sur un 
instrument : date d’achat, spécificattofîs,BédiiBÊ^e, réglage, réparations... 
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Exemple : Critères de vérification d’un spectrophotomètre 


• Alignement optique 

Vérification du centrage correct du faisceau lumineux ainsi que de sa stabilité 
(forme et taille) dans le temps : maintenance en atelier 

• Présence de lumière parasite 

Vérification de l’unicité de la longueur d’onde du monochromateur arrivant sur la 
cellule photoélectrique : maintenance en atelier 

• Exactitude photométrique 

Vérification de la justesse de la longueur d’onde affichée 

• Exactitude de l’absorbance et de sa répétabilité 

Vérification de l’absorbance connue d’un échantillon et de sa stabilité dans le 
temps. 


* Résolution spectrale 

Contrôle de la capacité du spectrophotomètre à séparer deux longueurs d’onde très 
proches (valable pour des spectrophotomètres haut de gamme) correspondant aux 
spécifications de l’appareil. 
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La métrologie d’un appareil est réalisée avec des 
matériaux de référence ou des étalons certifiés 


Etalons Primaires 

NIST 

N BS, N PL, LIME 




Certificats Officiels 
Nationaux 
ou Internationaux. 


Etalons Tertiaires 

Constructeurs 



1.2/ Validation de méthodes alternatives 


La validation d’une méthode alternative est basée sur : 

^ Méthodes de dosage dites de « référence » : 

• reconnues par la communauté scientifiques et faisant l’objet de NORMES. 

Méthodes réputées fidèles et justes et d’incertitude connue 
•mais souvent longues, chères, nécessitant des appareils sophistiqués... 


^ Utilisation par le laboratoire dans ce cas des méthodes « alternatives »: 

• utilisée à la place d’une méthode de référence. 

• pas forcément reconnue par la communauté scientifique ( pas de norme) 

• pouvant être propre à un laboratoire. 


Mais elle doit être VALIDEE, c’est à dire contrôlée. 
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Exemples de critères de validation de méthodes alternatives 

> Fidélité 


Une méthode est fidèle lorsqu’elle donne toujours le même résultat ou 
des résultats voisins si on la répète sur le même échantillon.. 


-Répétabilité: Fidélité sous des conditions de répétabilité (même méthode .même 
laboratoire, même opérateur, même équipement et pendant un court intervalle de temps) 


-Reproductibilité: Fidélité sous des conditions de reproductibilité (même méthode 

dans différents laboratoires, avec différents opérateurs et utilisant des équipements 
différents). 


> Justesse 


Etroitesse d'accord entre le résultat d’une mesure et la valeur attendue 
(CIBLE ou valeur réputée vraie) 
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Le centre de la cible 
représente la valeur 


attendue de la mesure 





Les points représentent les résultats de différentes mesures 
(en conditions de répétabilité ou de reproductibilité) 



Non fidèle 
Non juste 


Fidèle Non fidèle 

Non juste Juste 


Fidèle 

Juste 


Représentation de la fidélité et de la justesse sous forme de cible 
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> Sensibilité 


Rapport entre la variation instrumentale et la variation de la concentration. 


•Par exemple une méthode de dosage colorimétrique est sensible si elle est 
capable de donner des A franchement différentes pour des concentrations 
assez proches. 

> Limite de détection 

Plus petite quantité ou concentration distinguable d’un blanc. 


> Limite de quantification 

Plus petite quantité ou concentration pouvant être mesurée avec un risque d’erreur 
connu. 
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Les concepts de base de la métrologie en chimie 


la métrologie en chimie représente la voie à suivre pour améliorer la 
qualité des données analytiques. Elle englobe des concepts de base 
telles que : la traçabilité, l’incertitude, l’étalonnage et la validation des 
mesures, qui sont définies ci-dessous : 


A /Tracabilité 

Afin d'améliorer la qualité de l’analyse chimique; la métrologie est 
construite sur la possibilité d’établir une relation de comparaison entre 
la grandeur à mesurer, le mesurande et une référence connue de même 

nature. 

- Donc la traçabilité correspond à la possibilité de relier le résultat d’une 
mesure (ou de la valeur d’un étalon) à des références déterminées, 
habituellement des étalons internationaux ou nationaux afin d'obtenir 
un résultat quantifié. Pr r. bchitou 19 



B /L 'étalonnage : est l'opération qui permet d’effectuer des mesures 
de grandeurs connues Gi avec l'instrument de mesure donnant les 
valeurs Mi. On établit alors une courbe donnant les écarts entre les 
valeurs données par l'appareil et les valeurs des grandeurs connues. 

- On distingue deux types d’étalonnage: 

^ l’étalonnage de l’appareillage : qui correspond au test du paramètre 

physique mesuré par l’appareil. 

l’étalonnage de la méthode : qui permet d’établir la relation signal- 

quantité de substance ; ce dernier doit 
être réalisé avec des étalons et des 
matériaux de références adaptés. 


CI L Incertitude 

La valeur vraie d’une mesure ne peut jamais être connue, elle peut 
seulement être approchée. 

> L’incertitude sur la mesure est un paramètre associé au résultat de 
cette mesure, qui caractérise la dispersion des valeurs qui pourraient 
raisonnablement être attribuées au mesurage. 


D / Validation (Tune méthode d'analyse : est l’opération par 
laquelle on s’assure que les résultats répondent au problème de 
manière satisfaisante pour l’utilisateur. 

> Valider une méthode consiste à démontrer, avec un degré de 
confiance élevé et sous une forme documentée, que la méthode 
permet d’obtenir un résultat analytique qui atteint les spécifications 
définies à l’avance. 

>La validation se poursuit par une évaluation des résultats obtenus 
par la méthode pour l’analyse d’un matériau de référence certifié, de 
manière interne au laboratoire. 
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Système Internationale d’Unité (S.I) 


• Les unités SI sont donc, dans un système métrologique idéal, les 
références ou étalons primaires auxquels toute mesure doit pouvoir 
être traçable. 

•La métrologie est basée sur une référence internationale qui est le 
système international d’unités (SI) qui est composé des unités 
suivantes : 

- le mètre (m), unité de longueur ; 

- le kilogramme (kg), unité de masse ; 

- la seconde (s), unité de temps ; 

- l’ampère (A), unité d’intensité de courant ; 

- le kelvin (k), unité de température thermodynamique; 
- la mole (mol), unité de quantité de matière; 

Pr R. BCHITOU 22 

- la candela (cd), unité d’intensité lumineuse. 








• En chimie, le mole et le kilogramme sont les références 
primaires. 


•Les normes et les guides internationaux existants ont été 
établis selon des principes généraux et ne font pas de 
différences marquées entre les mesures physiques, 
chimiques ou biologiques. 
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Liste des organisations concernées par la métrologie en chimie 

ISO : International Organisation for Standardization. Genève 
BIPM : Bureau International des Poids et Mesures. 

CIPM : Comité International des Poids et Mesures. 

CCQM : Comité Consultatif pour la Quantité de Matière. 

CGPM : Conférence Générale des Poids et Mesures 

CITAC : Coopération on International Traceability in Analytical Chemistry 
SMT : Standard, Measurements and Testing Programme, Commission Européenne. 
Bruxelles 

REMCO : ISO Committee for Reference Materials. Genève 
NIST : National Institute of Standards and Technology. USA 

COMAR : International Data Bank on Reference Materials, Laboratoire National 
d’Essais. Paris 

CEN : Comité Européen de Normalisation. Genève 
AFNOR : Association Française de Normalisation. Paris 
ILAC : Internationl Laboratory Accréditation Conférence. 

EURACHEM : European Focus for Analytical Chemistry, Teddington. U. K. 

ICH : International Conférence on Harmonization of Technical Requirements for the 
Registration of Pharmaceuticals for Human Use. 

IUPAC : International Union of Pure and Applied Chemistry 

LNE : Laboratoire National d’EssaiSpP^ig CH|T0U 24 

FDA : Food and Drug Administration, USA. 


Exemple : 


Les normes se répandent dans l'industrie sont celles de la série des 
ISO 9000. 

les normes qualité concernant les laboratoires : norme ISO 17025 et BPL 


Des normes environnementales sont en cours de publication 
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Normalisation des instruments de mesure 

au Maroc 

Un comité de normalisation des instruments de mesure est 
institué afin de répondre aux besoins des industriels dans le 
domaine de la métrologie. Ce comité adopte des normes 
dans les domaines suivants: 

. généralités et définitions de la métrologie marocaine 
. mesure dimensionnelle 
. mesure de masse 
. mesure physico-chimique 
. mesure de température 
. mesure de pression 
. mesure de vitesse 

. mesure de volume et deRdéUtbu 2 e 




Chapitre II 



Pr R. BCHITOU 


27 




- 


' 

Comment présenter les résultats de mesures ? 


Encadrement des résultats 
(intervalles de confiance) 


Validité des résultats (tests d’hypothèses) 


Etablissement de prévisions (utilisation de 

modèles prévisionnels) 
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11.1/ Précision et exactitude 


Lorsque l’on réalise des mesures (des analyses), on veut 
rendre un résultat en s’interrogeant sur la validité de ce 
que l’on présente. 

(Plus la précision est grande , plus les indications sont proches de 
la valeur vraie). 

Dans ce contexte, cela implique l’évaluation des 
caractéristiques suivantes : 


► 

► 

► 

19 


Justesse (ou Exactitude) 

Fidélité (ou Précision) 
Justesse et Fidélité 
Robustesse ► Rugosité 



► Répétabilité 

► Reproductibilité 
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11.2/ Justesse, Fidélité et Erreur 

11.2.1/ Justesse ou Exactitude : Une méthode est réputée juste quand 
la moyenne x d'un grand nombre de mesures Xj est confondue avec 
la valeur X du mesurande , quelle que soit la dispersion . 

^L’erreur de justesse J est définie par : 

n 

ÜVeC X = X Xj/n ( moyenne des Xj) 

i=l 


J = X -X 


11.2.2/ Fidélité ou Précision : caractérise la dispersion d’une série de 
mesures X; d'une même grandeur (échantillon). Il est définit par l'écart 
type a ou la variance a 2 : 
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Fidélité (ou Précision) 



Pour un PROCEDE : 

à partir d’un échantillonnage 
du(des) produit(s) fabriqué(s) 
dans les conditions prescrites. 


Pour une METHODE : 

à partir de prélèvements multiples 
d’un échantillon homogène avec les 
conditions d’analyse prescrites. 


La fidélité doit être étudiée en utilisant des étalons ou des 
échantillons authentiques homogènes. 
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La fidélité peut être considérée à trois niveaux : 
Répétabilité, Reproductibilité et Précision intermédiaire . 

E Répétabilité 

( même série d ’ analyses ) 

l Reproductibilité 

(opérateur et jour et appareillage dijférents ) 

► Précision intermédiaire 

(opérateur, ou jour, ou appareillage différent ) 
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11.2.2.1/ Répétabilité 


Définition 


La répétabilité exprime la Fidélité pour les mêmes 
conditions opératoires pour une même série d’analyses 
dans un court intervalle de temps. 


V 

Elle est aussi appelée « précision intra-essai ». 
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II. 2. 2. 2 ! Reproductibilité 

Définition 


La Reproductibilité représente les Variations INTER-Ateliers 
ou INTER-laboratoires 


T 


C’est la mesure de la dispersion obtenue par plusieurs 
opérateurs qui analysent ou mesurent : 





C> dans des intervalles de temps différents, 


et 


> 


éventuellement avec des types d’appareils différents. 
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IL 2. 2. 3 /Précision intermédiaire 




la Précision intermédiaire, représente les Variations INTER- 
Ateliers ou INTER-laboratoires 



C’est la mesure de la dispersion obtenue par plusieurs 
opérateurs qui analysent ou mesurent : 



ou 


î> dans des intervalles de temps différents, 


ou 


■= > 


avec des types d’appareilsodifférents 
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11.2.3/ Les Erreurs Expérimentales 




elle varie toujours dans le même 
sens par rapport à la moyenne. 

i.e. une erreur qui reste constante 
pour des mesures effectuées dans 
des conditions identiques. 



Les erreurs systématiques 
affectent l'exactitude 


(Justesse) 


Pr R. 





elle se répartit de part et 
d'autre de la valeur moyenne 
(variance, écart-type ). 

J 



Les erreurs aléatoires sont 
relatives à la précision 
(Fidélité) 
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Recherche des causes d'erreurs 


Milieu Moyen Matière 



Mesure 


"les 5 M" 
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Démarche de Mesurage 



• Soit une série de n mesurages répétés. 



X 



3 » ■ 





, ,, , , ... Pi. H. BCHIT.OU , , 39... 

Idéalement, n tend vers I infini. Pratiquement, nous avons un échantillon 




Valeur 



vraie 
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11.3/ Intervalle de confiance 


Pour encadrer un résultat on parlera d’intervalle de 
confiance : 


0,9973 



Probabilité = 99,73 % pour que x 
soit compris dans l’intervalle p ± 3 a 



il y a plus de 99% de chances d’obtenir un résultat 
dont la valeur est égale à la valeur centrale +3o. 
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Exemple : 


Quatre opérateurs A, B, C et D dosent 10 ml de solution 0,1M de soude, mesurés 
exactement (précision instrumentale de 0,05 ml) avec une solution d'acide qui titre 
exactement 0,1 M . Quel est l’opérateur qui a travaillé avec précision et exactitude? 


A 

B 

C 

D 

10,08 

10,11 

10.09 

10.10 

10,12 

9,88 

10,14 

10,02 

9,80 

10,21 

10,19 

9,79 

9,69 

10,05 

9,78 

10,04 

9,98 

10,02 

9,97 

10,04 



moyenne 

10,10 

écart-type 

0,016 


Inexact^ Exact igg Inexact gg, Exact 


10,01 


0,172 


9,90 

10,01 

0,21 

0,033 


Précis 



Précis 
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Représentation graphique 



Valeur théorique de 10 ml, si on estime l’écart-type 
expérimental à 0,016 ml (opérateur A) alors : 

99,7 % des valeurs expérimentales doivent être 
comprises dans l’intervalle 10 ml + 0,048 ml 

Pr R. BCHITOU 
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11.4/ Méthodologie 


Pour la Maîtrise d’un Procédé ou d’une Méthode 
(d’analyse par exemple) il faut : 


•Procurer une connaissance totale et non biaisée des 
possibilités du Procédé ou de la Méthode telles que : 
justesse, fidélité et robustesse. 


•Structurer le travail expérimental de telle manière que 
les validations appropriées des caractéristiques du 
Procédé ou de la méthode puissent être considérées 
simultanément. 
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Procurer une connaissance totale et non biaisée des 
possibilités du Procédé ou de la Méthode telles que : 
justesse, fidélité et robustesse. 


I 


* 

I 






Structurer le travail expérimental de telle manière que 
les validations appropriées des caractéristiques du 
Procédé ou de la méthode puissent être considérées 
simultanément. 

Méthodologie des Plans d’Expériences 
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Voir E2 PIDD ET LIC PROF 
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11.4/ Mise au point de méthodes 

11.4.1/ ROBUSTESSE 


Définition 



La robustesse d’un procédé ou d’une méthode est une 
mesure de son aptitude à ne pas être affectée par de petites 
variations délibérées des paramètres de la méthode. 


Elle fournit une indication de sa Habilité pour un usage 
normal. 
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11.4.1.2/ Développement de Méthodes et Robustesse 


"L'analyse", "Le procédé », qu’ils soient chimiques ou physico- 
chimiques, impliquent la mise au point et l'utilisation de "méthodes' . 


Il peut s’agir : 


d'adapter une méthode existante au matériel dont on 
dispose ou à un nouveau type d'échantillons que l'on doit 
traiter (ajuster des volumes de réactifs, des temps et des 
températures de réaction et/ou des réglages d'appareils...) 
pour obtenir des performances satisfaisantes. 


^ de mettre au point une méthode originale 



se terminer p& f hlfê IT <8jkimisation 
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11.4.1.3/ Notion additionnelle à la Robustesse 


Pour la capacité d’un Procédé ou d’une Méthode à 
fournir des «produits» conformes on peut distinguer : 


La Robustesse 

4 * 

La Rugosité 

Faible sensibilité à une 

légère variation des 
facteurs expérimentaux 

maîtrisables. 

cr\ i ■ 

Faible sensibilité à une 

légère variation des 

facteurs expérimentaux 

non maîtrisables 


Paramètres du Procédé : Paramètres hors Procédé : 

température, concentration, temps, opérateur, espace, 
vitesse d’outils ... matériel, consommables... 
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II. 5/ Loi normale (ou loi de Gauss) 


Laplace et Gauss ont démontré que, pour la plupart 
des phénomènes physiques observables, les mesures 
expérimentales suivent une même loi de probabilité : 


une même fonction de densité de probabilité appelée 

Loi Normale. 


Pr R. BCHITOU 


51 




11.5.1/ Forme analytique de la loi normale 


Cette loi, qui décrit une variable aléatoire, 
] est caractérisée par deux 


Sa forme analytique est : 


un paramètre de dispersion : 
l’écart-type a. 


Un paramètre de position ou 
de centrage : la moyenne p 


paramètres : 




1 

t 

=L 

■ 

X 

v 

~2 

t 

Q 

V 
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II.5. 2/ Graphe de la Loi normale 


Distribution symétrique centrée sur la moyenne 


Point d’inflexion 
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11.5.3/ Propriétés de la loi Normale 



(fonction) = probabilité pour que la valeur de la variable X 

soit comprise entre -oo et x 


I 

Probabilité pour qu’une valeur d’abscisse soit 
comprise entre deux valeurs données ? 
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Probabilité p x pour qu’une 
valeur de x soit inférieure à x t 



Probabilité p 2 pour qu’une 
valeur de x soit inférieure à x 2 




p 2 - Pi = Probabilité pour qu’une 
valeur de x soit comprise entre x x et x 2 
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0,6827 0,9545 




Probabilité = 68,27 % pour que x soit Probabilité = 95,45 % pour que x soit 
compris dans l’intervalle p ± 1 g compris dans l’intervalle p ± 2 a 


0,9973 



Probabilité = 99,73 % pour que x 
soit compris dans l’intervalle p ± 3 a 
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0,6827 0,9545 




Probabilité = 68,27 % pour que x soit Probabilité = 95,45 % pour que x soit 
compris dans l’intervalle p ± 1 g compris dans l’intervalle p ± 2 g 


0,9973 



Probabilité = 99,73 % pour que x 
soit compris dans l’intervalle p ± 3 g 
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11.5.3/ Loi normale Standard 


Telle que nous venons de la définir, la loi Normale est 
fonction de p et a exprimés avec l’unité de la variable X : 


chaque cas est donc un cas particulier 



On peut rendre la loi universelle à l’aide d’un changement 
de variable : 



en prenant la moyenne p de la 
distribution pour origine de l’axe des x, 



avec l’écart type de la distribution 
comme unité de mesure. 
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Cette nouvelle variable s’appelle variable centrée réduite z, 

elle est sans dimension : 



Les caractéristiques de Z sont : moyenne = 0 et écart type = 1 


La forme analytique de la 
Loi Normale Standard est : 
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Loi normale standard 



C’est une loi universelle, indépendante des unités de la 

variable étudiée 


Elle s’utilise de la même manière que la loi normale 
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Intervalles de confiance 


Valeur individuelle : 

X ± Z G 


Grands Echantillons 
(> 30 répétitions) : 

X± z c ?- 

c \/n 


Petits Echantillons 
(< 30 répétitions) : 

x±t c ® 

v n 
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11.5.3/ La distribution de Gauss ou « normale » 


Si l'on réalise un nombre infini de mesures sur le même échantillon et dans les 
mêmes conditions expérimentales, on obtient une série de résultats qui se 
répartissent avec des fréquences qui suit une distribution normale. 


série de n mesures : y ,\ x 2 ... 

I xi 

moyenne |i = = médiane 

n 


2 1 » . 2 
variance = G = ” L ( xi- ju) 


n . 


= 1 


écartype = G 




95 , 5 % 
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moyenne différentes 
dispersion identique 



•Plus l'écart-type est grand, plus la dispersion est grande 


Répartition des mesures, intervalle de « confiance » 


A 

éf 


a 


/ 

/ 

/ 




r 

+ 

-e- 

n 

la 

X 


environ 68% des mesures 
sont comprises dans 1 intervalle 
m ± la 




environ 99,8% mesures 
sont comprises dans 1 intervalle 
m±3a 63 


II. 6/ Modèle statistique des essais inter-laboratoire 

r'Le laboratoire L, effectue un dosage sur l’échantillon E avec une méthode donnée 
en condition de répétabilité , il obtient une série de mesure qui suit une distribution 
gaussienne: moyenne m b écart type de répétabilité a rl 




^Le laboratoire effectue un dosage sur le même échantillon E avec la même 
méthode et en condition de répétabilité, il obtient une série de mesure qui suit une 
distribution gaussienne : moyenne m i5 écart type de répétabilité a ri 


Si la méthode est bien normalisée, les écart types de répétabilité a ri sont 

identiques. C’est cette valeur identique qui est nommée l’écart type de 
H Pr R. BCHITOU 
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64 



* Mais les moyennes des laboratoires sont différentes 




valeur de x 


Si on fait la moyenne des moyennes : m= (mi+m 2 + .m n )/n 


Le jour J, un laboratoire donné L; donne un résultat qui s’écarte de 
la moyenne m avec une déviation B ( biais du laboratoire le jour J) 


B prend une valeur aléatoire, vafte R tëüB l Té& J jours, sa moyenne est nullé 5 


- En reportant les moyenne de chaque laboratoire sur un histogramme, on obtient une 
distribution normale de moyenne m et d’écart type s interlaboratoire 



•Soit X le résultat d’un dosage, effectué par un laboratoire L i5 un jour J 


Valeur supposée vraie 
Mathématiquement : 


x = m + B + e 


Biais de ^ le jour J 




Aléa de répétabilité de la 
méthode 

y 


Aléa de reproductibilité 



^ R ^ reproductibilité — \/ (c>,nterlaboratoire) + (^répétabilité) 


D’où 


Variance ( B+e) ^r\ferèaiiMi:eo(iB)+ Variance (e) 66 

OU (o reproductibilité) ” (^interlaboratoire) + (^répétabilité) 



M 

Mais quelle est la valeur « vraie » du dosage? 

^La valeur vraie n’est connue que si l’on dispose d’un matériau de référence ( MR). 
C’est à dire un échantillon dont la valeur de la mesure est parfaitement connue 

>^Les matériaux de référence peuvent être certifiés ( MRC). C’est à dire accompagné d’un 
certificat garantissant la valeur réelle: jn. Ce sont des étalons garantis. 

>^La méthode de dosage peut alors donner un résultat différent de p, avec une 
déviation 8 ( biais de la méthode). 

Si on ne dispose de MRC, on admet que la moyenne des moyennes donnent la 
valeur vraie: m est confondue ave^jLi R bchitou 67 


Remarque importante 


•Avant tout calcul statistique, il est nécessaire d’éliminer les valeurs 
aberrantes, sinon la moyenne, l’écart type n’ont plus de sens... 

•Pour éliminer les aberrants, on réalise des tests mathématiques ( que 

nous n’aborderons pas dans ce cours). 

Test du x 2 : pour comparer les valeurs x, intra-laboratoire et rejeter les 
aberrants avant de calculer moyenne et écart type. 

Test de Cochran: pour comparer les écart types s ri et rejeter les 
aberrants avant de calculer s r 

Test de Grubbs: pour comparer les moyennes ïw { et rejeter les aberrants 
avant de calculer m 

Observez les résultats avant de vous lancer dans les calculs 


Métrologie en chimie 


Quelqilfé^pWMpëé et beaucoup d’incertitudes 


Chapitre III 
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LES PLANS D’EXPERIENCES 


Outils indispensables pour le développement des Méthodes 
analytiques : 

- au niveau du procédé de préparation de 
l’échantillon analytique 

- au niveau de l’optimisation de la méthode 


D’une manière générale : à mettre en œuvre au niveau de 
toutes les étapes expérimentales de développement et de 
validation de la méthode. 
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Expériences 




Phénomène 
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Plans d 'expérience 


Facteurs 


Réponses 





Modèle 

mathématique 






Yl 

y 2 


variation n 

Raisonnée Calcul des Effets 

des Facteurs sur les réponses 


Pr R. BCHITOU 


72 






Rappel de quelques définitions concernant 

Les plans d'expériences 


Matrice d'expérience : 

Matrice d'expérience est un tableau de n lignes et k 
colonnes, regroupant les conditions expérimentales d'un plan 
d'expériences, n et k correspondent respectivement au nombre 
d'expériences et au nombre des variables codées. 

Variables explicatives et notion d’interaction : 

Les variables explicatives d'une étude sont les paramètres 
susceptibles de modifier les réponses de cette étude. Si l'effet 
d'une variable explicative dépend du niveau d'une autre 
variable explicative, on dit qu'il y a interaction entre ces deux 
variables explicatives. 
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Niveaux d'une variable explicative : 

Les niveaux d'une variable explicative sont les différents 
états que peut prendre cette variable explicative 

Notion d'effet significatif : 

L'effet d'une variable explicative sur la réponse y 
s'obtient en comparant les deux résultats de mesure y, et 
y 2 de réponse, mesurée lorsque la variable explicative 
passe d'un niveau (0) à un niveau (+). Si l'écart entre y t et 
y 2 est important on dit que le facteur est influent ou 
significatif. 

Variables codées et variables naturelles : 

Les variables naturelles xi sont les valeurs qui 
correspondent à chaque niveau d'une variable explicative. 
Pour comparer les effets des variables naturelles sur la 
réponse, il est nécessaire de les remplacer par les 
variables codées Xi qui sont sans unité. 
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Réponse : 

Le résultat mesuré d'une étude. A chaque point du 
domaine d'étude correspond une réponse. 
L'ensemble des réponses forme la surface de 
réponse. 

Courbes d'isoréponses : 

Après la détermination du modèle et la 
vérification de sa validité, les courbes d'isoréponses 
peuvent être tracé à l'intérieur du domaine 
expérimental. Ces courbes représentent des plans 
pour surfaces de réponse c'est à dire la 
représentation graphique des résultats (modèle 
estimé) pour pouvoir en tirer des optimums. 
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Variance-covariance 



Variance de x 




1 


n 


°y = 


n-1 w 


Z(K-7) 


Variance de y 




Covariance x 1 
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111.1/ MATRICES D’ EXPERIENCES et 
METHODOLOGIE 

Les objectifs : 

• Le criblage des facteurs : classement hiérarchisé des 
facteurs; 

• Les études quantitatives des facteurs : quantification des 
influences principales et des synergies éventuelles; 

• Les études quantitatives des réponses : modélisation 
prévisionnelle du phénomène étudié; 

• L'optimisation : déterminer un ou plusieurs points de 

fonctionnement optimaux . 
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III .1.1/ Les Modèles : polynômes 


Soit X un facteur quantitatif. Il peut être représenté par un polynôme 
dans les modèles : 



Premier degré y (réponse) = b 

>D + b , *1 + b 2 X 2 


i— ^ / 

Effet principal! 


Synergique 


Y (réponse) = b 0 + b., X 1 + b 2 X 2 + b 12 X,X 2 


\ 


Effet d'interaction 


Quadratique 
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III .1.2/ Les Matrices 


MATRICES 


Premier degré : 


Y (réponse) = b D , + b, X 1 + b 3 X 3 


^=> 


Hadamard 

(Plackett-Burman) 



Y (réponse) = b D + b, X, + b 2 X 2 + b 12 X^ 



Factorielles 


Quadratique : 


Y (réponse) = b a + b 1 X! + b 2 X 2 + b 13 XJ( ; + b u X 2 + b 2a X 2 2 


Composite 

(Doehlert) 
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III. 2/ MODELISATION 


Modéliser : utiliser des données expérimentales pour prévoir une 
information quantitative inconnue Y à partir de mesures de X via une 
certaine « fonction mathématique » : 



Le modèle mathématique postulé peut être : 




Ajustement polynomial 


Moûm Y = ftX)e.*c un 
povnôme se Segre 
MMMMl * 


Une droite si Y varie 
linéairement avec X. 
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Exemple : Etude de la stabilité d’une suspension. 

> Chercher un modèle liant les facteurs expérimentaux. 


Facteurs expérimentaux 

Variables 

Domaine Expérimental 

Minimum 

Maximum 

Tensio-actif : Mouillant 1 

M, 

20 g/l 

40 g/l 

Tensio-actif : Mouillant 2 

m 2 

5 " 

15 " 

Epaississant : Structurant 1 

s, 

5 M 

20 " 

Epaississant : Structurant 2 

Sj 

0 *' 

10 " 


Réponse étudiée : 

Y- c o de séparation de la suspension en deux phases 
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Tableau des résultats 


Plan d'expérimentation 


Réponse 


Essais 

Mi 

m 2 

Si 

s 2 

Y 

1 

40 

15 

20 

0 

10 

2 

40 

15 

5 

10 

16 

3 

40 

5 

20 

0 

15 

4 

20 

15 

5 

0 

38,7 

5 

40 

5 

5 

10 

30,5 

6 

20 

5 

20 

10 

18 

7 

20 

15 

20 

10 

13 

8 

20 

5 

5 

0' 

32 


Modèle postulé 


Y = a c , + a 1 + a 2 a 4 S 2 
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III. 2. 1/ Codage des variables : 

-Variables naturelles : U 
-Variables codées : x 


Ui 

Toutes les variables 
naturelles ont leurs propres 
plage de variation qui 
dépendent des unités 

Les coefficients du modèle ne 
peuvent pas se comparer 
directement 

Interprétation difficile 


Xi 

Toutes les variables codées 
ont la même plage de 
variation entre -1 et +1 
indépendante des unités 

Les coefficients du modèle 
peuvent se comparer 
directement 

Interprétation facile 
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= moyenne des valeurs maximum et 
minimum que peut prendre la variable Uj 


U:: = valeur de la variable naturelle 1 à l'expérience i. 


ij 





u„ ■ u.° 


AU, 


U 1 -, = valeur de lia variable naturelle j 


au centre du domaine. 



AUj = pas üe variation 
de là variable naturelle j. 


X;, = valeur de la variable codée j pour l'expérience i. 


= demi étendue : moitié de l'écart entre la 
valeur maximum et la valeur minimum que 
peut prendre la variable Uj . 


Réciproquement: 

U, = IV ♦ X, . AU, 
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Valeurs codées 


est le mouillant 1 qui varie entre 
1T 20 et 40 a /l ^ 

centre du domaine l). 

Pas de variation AU, 

= (20+40}/2 = 30g:l 

= (40 -20) *2 = 10 g/l 

ÇU à 20 g/l: X, t = {20 - 30)/1 0 
^ à 40 g.'l: X M = (40 - 30>/10 

ü ü 




U4 es! le structurant 2 qui varie entre 

0 et 1 0 g.'l 


centre du domaine U a 4 


Pas de variation MJ, 

= {0+10J/2 = 5 g;l 


= (1 0-ÜJ/2 = S g/l 

à 0 g/l : X M 
àlO g/l : X u 

= (0-5)/5 

= -i 

,= {10-5J/5 

= i ^ 
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Régression avec des variables codées 


Matrice d’expériences 


Essais 

Xi 

X 2 

Xi 

X, 

1 

+1 

-I 

+1 

-1 

2 

+1 

-I 

-1 

+1 

3 

+1 

-I 

+1 

-1 

4 

-1 

-I 

1 

-I 

5 

+1 

-I 

-1 

+1 

fi 

-1 

-1 

+1 

+1 

7 

-1 

+1 

+1 

+1 

S 

-1 

-I 

-1 

-I 



Modèle postulé 

Y = b c , + bi Xi + b 2 X 2 + b ? X 2 + b 4 X 4 


Comment déterm4i$f 0 î>o> b i? b 2 > b 3 et b 4 ? 86 


Chapitre IV 
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Balance 



Nous convenons d’appeler ••expérience” le fait cle lire la position 
de l’aiguille sur le cadran. 


Comme pour toute expérience, il existe une erreur expérimentale : 


erreur de parallaxe, épaisseurs relatives de la pointe de 
l'aiguille et des traitspideEaRtdMation etc. 88 


Efficacité 


r 


Etre efficace c'est rendre 
minimum i 



non pas Teneur de lecture 


MAIS 


sa répercussion sur 
l'estimation des masses, 


l 

nous admettrons que 
i'expé ri me ntate ur pre nd! to utes 
les précaytions utiles pour la 
rendre minimale 


i 

(sur le calcul). En d’autres 
ternies c'est cerner au plus près 
la masse vraie des objets à 
peser. 
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Lecture avec la pomme 


Masse de la pomme = y 2 - y. 
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Chaque pesée, matérialisée par la 
lecture Y (réponse), est : 


indèpen dante des autres u n e va r iab I e a I éato i re car Y, est 

lectures une mesure expérimentale. 




Y, est en réalité la somme de deux valeurs : 

^i = Tii + e| 




n | = valeur vraie 

e, = erreur de mesure aléatoire 
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Caractéristiques de l'erreur expérimentale r 


Distribution de Gauss centrée sur zéro 
En moyenne, # 0 
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Qualité de l 'estimation de la masse de la pomme 


l'estimation de la Masse de la pomme m = y 2 - y, 

var(m) = var{Y 2 - Y.,} = var(Y 2 } + varfY,) 

= var{ri 2 + e 2 ) + varfr, + ej 
= var(Ti 2 )+ var(e 2 ) + var(ri ,) + var(e,) 

— n x / — 

“0 = CF 2 

(variance de e ) 

var(m) =2 o 2 
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Pesée de plusieurs objets 


Si nous avons trois objets, comment les peser ? 
Quel est le prix minimum de l’expérimentation? 
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Expérimentateur n°1 




Masse d’un objet : m- = yj., - y, 



var(fnj) = var(YjJ ♦ var(Y,) 


var(m i ) = 2 a 2 
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Matrice d ’expériences 


Une matrice d’expériences est: un tableau permettant de décrire 
une expérimentation en donnant, pour chaque expérience, les 
valeurs des facteurs expérimentaux. 



d ,j caractérise l’état de l'objet j 
au cours de l'expérience i 
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Expérimentateurs n°1 & n°2 


Codage des objets 


. objet absent : état 0 
' objet présent : étal 1 


Exp. 

0 

i Cl -jji 

0* 

1 

0 

0 

0 


2 

1 

1 

0 


3 

1 

0 

1 


4 

0 

1 

1 



Plan N°2 (D s ) 


Exp. 

0 

1 °2 °3 

1 

“o 

0 0 


2 

1 

0 0 


3 

0 

1 0 


4 

0 

0 1 



Plan N°1 {D,) 
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Expérimentateur n° 2 


TT 



Yi 




I 


iLectu res 


y 2 ) <*T 1 ' 




y a 


_ ï" 

a 


V\ 


\ 





Q 


<r r r > ^ 



j 


Système d'éq u alio us 


ïi 

y* 

y d 

y* 


= Hj + *i 

= n, + n* 
= n, + n a 
= ï: + n 3 


= m 5 


= 

e< - 


e, = 


- _r -p 

m ü + m 1 * rrr 
m u + in., + m, 
m., + itv + fn 3 
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Matrice d 'expériences et système d 'équations 
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Matrice du modèle et système d équations 


Y = m (l + nifl X 1 + m £ X 2 + m 3 X 3 

Avec la notation matriciel le :r oni peut écrire de manière compacte 
le système de 4 équations linéaires qui lient les réponses 
observées aux états des obj ets : 

Y = Xm 






IL 



us 


1 

-1 -1 1 


m D 






1 

1 -1 -1 

■* 



Vecteur réponse 

Vs 


1 

-1 1 -1 


m 2 




Vt 


1 

1 1 1 


m 3 












Matrice du modèlle 


Vecteur elles coefficients 
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Expérimentateur n°2 




Résolution 


*»o = >'i 


Lectures 



= (- > 1 + y 2 + y 3 - y 4 ) / 2 
in> = (-yi + y 2 -y 3 +y 4 )/2 

‘“3 = (->'i-y2 + >'i + >'4) /2 

var(mi) = var [l^X^. de 4 Yj)J 
= 1/4 [var(Z algibr de 4 Y,)] 





Expérimentateur n°2 




Résolution 


*»o = >'i 


Lectures 



= (- > 1 + y 2 + y 3 - y 4 ) / 2 
in> = (-yi + y 2 -y 3 +y 4 )/2 

‘“3 = (->'i-y2 + >'i + >'4) /2 

var(mi) = var [l^X^. de 4 Yj)J 
= 1/4 [var(Z algibr de 4 Y,)] 





Expérimentateur n°3 


□ 

/ Codage des objets : 

àsfc, J 

' 3 

i 

i 

objet présent à droite : état 1 

C 1 

il 

1 

> 

-4- 

Exp. O | 0 2 0 3 

i r-i -i -i 

£ 1 -1 -1 

3 -1 1 -1 

4 |_-1 _1 1 

Plan NT3 (D 3 ) 


4 . 
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Expérimentateur r 

i°3 




ÆJbi 



Système d'équations 


i - ■ irijj - iiij - uj. > - nij 

y 2 = fi 0 + m 1 - m\ - ïn J 
Yi = in 0 - fil + fi; - fij 


y 4 = m 0 - mi - m 2 + mj 
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Expérimentateur n°3 



Lectures 



Résolution 

m 0 = (-yi + y 2 + y 3 +y 4 ) /2 

“1 = (>'2 - >1) ' 2 

m 2 = (y 3 - >i) / 2 

û 3= 0*4 - >1) / 2 



var(mj) = vai[l/2(I al?tbr de 2 Yj)] 
= 1/4 [vard^ de 2 \\)] 

var(m i ) = a 2 /2 
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Expérimentateur n°4 

n 



\ 

^ Seule l'expérience 1 change 
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Expérimentateur n°4 





x 


Système d’équations 


Yj = id i;i + m l + hgl 2 + ni 3 
y 2 = - m 2 - 


Ya = mo - m i + m 2~ ms 


j j = m 0 - in L - in 2 + în^ 
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Expérimentateur n°4 


Lectures 



Résolution 

“o - (+ >1 + >'2 + > 3 +>'4 ) / 4 
“i“( + yi + y2-y*-y4)/‘* 
ô 2“( + yi-y2 + yi-y4> /4 
= ( + >i - y 2 - y 3 + y 4 ) / 4 

var(m|) = var[l/4(Z alsibr de 4 Yj)J 
= 1/16 [vai e:^^ de 4 Y,)) 



varCnii) = a 2 /4 


L 


X 


J 
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Conclusion 
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THÉORÈME 

Si on fait N pesées, la variance sur les coefficients est > cf 2 'N 


Comment obtenir une matrice optimale 
pour minimiser la variance sur les masses estimées? 
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Modélisation 


La relation entre la mesure Y (la lecture du cadran) et les différentes 
masses dont on cherche à déterminer la valeur est une relation 
linéaire qui s’écrit : 


masses des objets 
m 0 , m 1f m 2 et m } 



état des masses dans erreur 
l'expérience j (la pesée) 

X 1 . X 2 et Xj présence ou absence 

^M^ ux (- 1 - 1ou0 »- 
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Matrices d ’Expériences 


Exp. 

o, 

O3 

1 

r° 

0 on 


2 

1 

0 0 


3 

0 

1 0 


4 

0 

0 1 



Plan N°1 (D,) 


Exp. 

0, 

o 2 

o 3 

1 

-1 

-1 

-1 

2 

1 

-1 

-1 

3 

-1 

1 

-1 

4 

-1 

-1 

1 


Plan N°3 (D 3 ) 


Exp. 

o, 

o 2 

°3 

1 

0 

0 

0 


2 

1 

1 

0 


3 

1 

0 

1 


4 

0 

1 

1 



Plan N°2 (D 2 ) 


Exp. 

o. 

o 2 

o 3 

1 

1 

1 

1 


2 

1 

-1 

-1 


3 

-1 

1 

-1 


4 

_-1 

-1 

1 



Plan NM (D 4 ) 
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Expression matricielle de la régression 


Y = m 0 + m 1 X 1 + m 2 X 2 + m 3 X 3 + e 

Avec la notation matricielle, on peut écrire de manière 
compacte le système de 4 équations linéaires qui lient les 
réponses observées aux états des obj ets : 

Y = Xm + e 


# 

M 




y 

c J 

m 

L 


X 





Yi 


1-1-11 


m o 


e. 





v 2 


i i-i-i 

* 

m 1 


e 2 


Vecteur 

Vecteur réponse 

Ys 


i -i i -i 


m 3 

T 

e s 


des erreurs 


y* 









i n i i 


Jll 3 







— 










Matrice du modèle 


Vecteur des coefficients 
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Expression matricielle de la régression 

La matrice X représente lia matrice des coefficients des 
inconnues in, du système d'équations. 


Calcul du vecteur des coefficients 


X' étant la matrice transposée de X. si X’X n'est pas singulière 
{c'est à dire s’il est possible de calculer la matrice inverse 
correspondante) : 


on aura l’estimation du vecteur des coefficients du modèle au 
sens des moindres carrés par : 


m = [X , X)' 1 X , Y 
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Qualité des Matrices d Expériences 


EX|>™ 

1 

2 

3 

4 


Oi o 2 o 3 

' 0 0 G 1 

I 0 0 

0 I 0 

, 0 0 1 , 

Matrice D 



Matrice de variance-covariance : (X’X) _1 
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r 

Analyse quantitative et Etalonnage 



Les données sont toujours en nombre limité 


>Elles ne représentent donc qu’un échantillon de la population de 
toutes les mesures de la teneur en analyse de l’étalon que l’on 
pourrait effectuer. 



- Si X représente une teneur Y représente le résultat 

connue en analyse observé 


'~ J > La relation linéaire postulée devient : Y= b 0 +b 1 X 



Avec uniquement une "estimation" 
des coefficients b 0 et b x du modèle postulé. 

Ÿ= p 0 +PiX 


Pr R. BCHITOU 


117 









Expression matricielle 


( 


! 


1 

x 2 

1 

x 3 

1 

x 4 

k ■ ■ 

1 

1 

k è S 

i 


* 





Yi 

y 2 

y 3 

y 4 



Vecteur des 
résidus : T 


Vecteur de la 

réponse 

expérimentale 
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Estimation des coefficients : 


Les (3j sont les inconnus que nous devons estimer : 
(bj est l’estimation calculée de p, ) 


1 . Au sens des moindres carrés (résolution algébrique) : 
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2. Au sens des moindres carrées (résolution matricielle) : 


>En notant X' la matrice transposée de X, on aura l’estimation du 
vecteur des coefficients du modèle au sens des moindres carrées par 


Vecteur des coefficients 


(X'X> <X'Y 
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Chapitre V 
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But : 


- Le criblage des facteurs : classement 
hiérarchisé des facteurs 
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Effets des Facteurs 



Effet de X, : E, 


Ei = Y b -Y a 


Y 

Eip.V 


il 

1 

-1 



1 


Y 




Exp.: 



+1 


X 


Effet de X : 

(va eu r de a réponse Y au niveau + 1 ) 
iValeur ce la réponse Y au niveau -1 ) 


Pr R. BCHITOU 


123 


Poids des Facteurs 



Exp„V 3 


1 









d u < entre du domaine 
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Poids des Facteurs 



Exp.X 5 



1 


-1 

■) 


-1 

è 


1 

J 


1 




Modèle : Y = b„+ b, X,+ b„X, 


C jL 




de X, l > 1 


Ei.pt-> 

L 

K\p.4 

-1 

41 

1 








FYp l 


F.ipwï 

> i + 



y* 



Valeur moyenne (> 


N , + \ , + \ -.Hh 

b I fa fa -I fa 


il 


Poids = Pente = Coefficient 
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Modèle pour un facteur : Y = b 0 + b t X 
Modèle pour deux facteurs : Y = b 0 + b 1 X 1 + b 2 X 2 

Modèle pour K facteurs : Y = b 0 + b,X, + b 2 X 2 + b 3 x 3 ... + b k X k 


Pour avoir la variance des estimations minimale a 2 / N, 
la condition nécessaire et suffisante est que l'on ait : 


X'X = N l N 





Plans de Pesées et Matrices d ’Hadamard 


Cette propriété est obtenue avec les matrices 
d'Hadamard, qui ne contiennent que des 
valeurs -1 ou +1 . 

Elles n'existent que pour N multiple de 4 : 

4. 3, 12, 16, 20 ... 



Jacques Salomon 
HADÀIMIÀRD 
- 1963-} 
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Construction de Plackett et Burman 

k < 3 (N =4) : + + - 

k <7 (N=8) : + + + -# -- 

k < 11 (N=12) : + + - + + + — + - 

k < 15 (N=16) : + + + + - + - + + -- + — 

k < 19 {N=20) : + + + + + + - + + 

k < 23 (N=24) : + + + + + -+ -+ +- - + +-- + - 


Exemple N = 4 

® + - 
+ - 


+ + + - 

après N -2 permutations 

- + + 

■ "1 !■"■ 


+ - + + 
+ - - - 


Matrjç$fi c 

HITOU 

Matrice X 








Plans de Pesées et Matrices d ’Hadamard 


Transpositl 



+ 

+ + + 


+ 

+ - 


+ 

_ + 


- 

+ + 


Patrice X : 



4 0 0 0 

0 4 0 0 
0 0 4 0 

0 0 0 4 

Patrice (X'X) 


1/4 

0 

0 0 

0 

1/4 

0 0 

0 

0 

1/4 0 

0 

0 

0 1/4 

Patrice 

(X'X}- 1 
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Matrice d'expérience 


k<7<N=8) : + + + - + -- Décalage droite 


N* Esp 

XI 

X2 

X3 

X4 

X5 

X6 

X7 

1 

i 

i 

i 

4 

i 



2 


i 

i 

1 

-i 

i 

-i ( 

2 








4 








5 








S 








7 








2 
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Matrice d'expérience 


k < 7 (N=8) : + * + - + Décalage droite 


N= Exp 

Xli 

X2 

X3 

X4 

X5 

X6 

X7 

1 

i 

i 

i 

-i 

i 



2 

-i 

i 

-a 

i 

-i 



3 




i 



* l 

4 








5 








S 








7 








3 
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Matrice d’expérience 


k<7{N=8) : + + + - + -- 


N 3 Eïp 

XI 

X2 

X3 

X4 

X5 

X6 

X7 

1 

i 

i 

i 

-i 

i 

-i 

-i 

2 

-i 

i 

i 

i 

-i 

i 

-i 

3 

-i 

-i 

i 

i 

i 

-i 

i 

4 



-i 

i 

i 

1 

-i 

5 

-i 

i 

-i 

-i 

i 

1 

i 

G 

i 

-i 

i 

-i 

-i 

1 

i 

7 

i 

i 

-i 

i 

-i 

-1 

i 

3 

-1 

-1 

-i 

-1 

-1 

-1 

-1 


Dernière ligr 

ip 1 ________ 

IC . 

Matrice complète 
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Chapitre VI 


Plans factoriels complets 

FFD 
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Les plans factoriels complets sont caractérisés par : 



Pour k = 2, on obtient : 

> Une matrice d’expériences 2 2 : soit N= 4 lignes et 2 colonnes. 
r La matrice du modèle X : N= 4 lignes et 4 colonnes. 
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Modèle synergique : deux facteurs 


Plan factoriel complet pour 2 facteurs avec interaction 



X , 

x a 

-1 

-1 

1 

-1 

-1 

1 

1 

1 



Matrice d’expériences 
Plan factoriel complet 2 2 

Toutes les combinaisons de 2 facteurs 
pouvant prendre 2 niveaux 


|ModèlejY = bg+b^Xj+b^X^ 


] 







Exp.N 3 

x„ 

X, 

X, 

X,x 2 


1 

1 

"-1 

-1 

1 


2 

1 

1 

-1 

-1 


3 

1 

.1 

1 

-1 


4 

1 

1 

1 

1 



Matrice du modèle 
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Matrice d’expériences: 3 facteurs 


Essai 

XI 

X2 

X3 

1 

-1 

1 

-1 

2 

1 

1 

-1 

3 

-1 

1 

-1 

4 

1 

1 

-1 

5 

-1 

1 

1 

6 

1 

1 

1 

7 

-1 

1 

1 

B 

1 

1 

1 


Matrice d'expériences 

Plan factoriel complet 2 S 

Toutes les combinaisons 
de 3 facteurs pouvant 
prendre 2 niveaux 

Algorithme de Yates 


(X3) 

A 


FFD - 

effets / interactions 
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Modèle synergique : 3 facteurs 


Y = b 0 + b 1 X 1 + b 2 X 2 + b 3 X 3 + b 12 X 1 X 2 + b 13 X-,X 3 + b 23 X 2 X 3 

V A- *ÿl23 X 1 X 2^3 J_ 


/ 3 effets d’interaction 

1 effet d’interaction du premier ordre b 
du second ordre b ijk 


1 terme 3 effets 
constant principaux bj 


N ° exp . 

Xo 

X, x 2 

X 3 

x,x 2 x,x. 

x* 

X,x 2 x 3 

1 

1 

-1 -1 

-1 

1 1 

1 

.1 

2 

1 

1 -1 

-1 

-1 -1 

1 

1 

3 

1 

-1 1 

-1 

-1 1 

-1 

1 

4 

1 

1 1 

-1 

1 -1 

-1 

-1 

5 

1 

-1 -1 

1 

1 -1 

-1 

1 

6 

1 

1 -1 

1 

-1 1 

-1 

-1 

7 

1 

-1 1 

1 

-1 -1 

1 

-1 

8 

1 

1 1 

1 

1 1 

1 

1 
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Matrice du modèle 
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En conclusion, pour k=3, on obtient : 


Matrices d'expériences : matrice de taille 2 3 
soit N=8 lignes et 3 colonnes 

Matrice du modèle X : 8 
lignes et 8 colonnes 


N°exp. 

X 0 

X 1 x 2 x 3 

x,x 2 

X,x 3 

x 2 x 3 

X,X 2 X 3 

1 

1 

-1 -1 -1 

1 

1 

1 

-1 

2 

1 

1 -1 -1 

-1 

-1 

1 

1 

3 

1 

-1 1 -1 

- 1 

1 

-1 

1 

4 

1 

1 1 -1 

1 

-1 

-1 

- 1 

5 

1 

-1 -1 1 

1 

-1 

-1 

1 

6 

1 

1 -1 1 

- 1 

1 

-1 

-1 

7 

1 

-1 1 1 

-1 

-1 

1 

-1 

8 

1 

1 1 Pr R. 4 ch 
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1 
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Exemple d'utilisation dune matrice factorielle complète : 
étude quantitative des facteurs 

La réaction enzymatique est réalisée par la culture de la 

souche microbienne sélectionnée, dans un milieu nutritif 

contenant le substrat à déshydrogéner : 

c'est le développement de la souche qui libère le système 

enzymatique dans le milieu nutritif et qui "métabolise" le 

substrat. 


Facteurs 

Bornes 


(-D 

(D 

U 1 : Quantité de liqueur de maïs en g/l 

10 

20 

U 2 : Durée de réaction en heures 

24 

43 

U 3 : Quantité de glucose en g/l 

5 

10 


RÉPONSE Y étudiée : 


Quantité de substrat (mg) d^sfh^elfïQ^née enzymatiquement . 141 



Etude d T une déshydrogénation enzymatique 


Matrice d’expériences Plan d’expérimentation 



Ui 

Quantité 
.mais 


20 

10 

20 

10 

20 

10 

20 


U 2 

Durée 

réaction 


24 


48 

48 

24 

24 

48 

43 


Oté 

glucose 


5 

5 

5 

5 

10 

10 

10 

10 


mg 

subst 

transf. 


230 

205 

110 

70 

270 

220 

110 

70 
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Matrice du modèle 


Exp . N 0 

x 0 

Xi 

x 2 

X 3 

XiX 2 

XiX 3 

X 2 X 3 

XiX 2 X 3 

Y 

1 

1 

-1 

■1 

-1 

1 

1 

1 

-1 

230 

2 

1 

1 

-1 

■1 

-1 

-1 

1 

1 

205 

3 

1 

■1 

1 

-1 

■1 

1 

-1 

1 

110 

4 

1 

1 

1 

-1 

1 

■1 

-1 

■1 

70 

5 

1 

■1 

-1 

1 

1 

■1 

-1 

1 

270 

6 

1 

1 

■1 

1 

-1 

1 

-1 

■1 

220 

7 

1 

■1 

1 

1 

-1 

-1 

1 

■1 

110 

8 

_1 

1 

1 

1 

1 

1 

1 

1 

70 



l 2 

^3 

«4 

^5 

h 

II 

h 
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Exploitation des résultats : 


• a / Calcul des coefficients ( Déterminer une 
estimation ponctuelle des effets principaux sans et 
avec interaction ). 

• b / Analyse Statistique des résultats 

• c/ Significativité des coefficients 
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Matrice d’expériences 


X, 


-1 

1 

-1 

1 

-1 

1 


x 2 

X, 

Y 

T™ 

“1 


205 v 

~r~ 

-1 

-1 


m 

70 


■1 


1 


270 

220 


1 

n 

i 

110 

1 

i 

i 

^ m 

70 
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a / Calcul de coefficients 


Nom 

Coefficients 

b 0 

160,6 

bi 

-19,4 

b 2 

-70,6 

b 3 

6,9 

b 12 

-0,6 

b 13 

-3,1 

^23 

-6,9 

b 1 23 

3,1 


La durée de lia réaction (facteur X 2 ) a 
le plus d'influence : en moyenne le 
rendement diminue de 140 g/l quand 
on passe de 24 à 43 h. 

Cela signifie que le produit de 
déshydrogénation est à son tour 
métabolisé, il ne faut pas laisser la 
culture trop se développer. 

La quantité de liqueur de mais 
(facteur XJ doit également être 
choisi à son niveau minimum (10 g/l), 

La quantité de glucose (facteur XJ 
n'a proportionnellement que peu 
d'effet. 
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c/ Significativité des coefficients 



b i" b i° 
t= é.t.(b) 

t=-^- 

é.t.(bi) 
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Exploitation des résultats 

Réexamen «statistique» en négligeant b 123 comme coefficient : 

r — v JL 

j /sa valeur est alors considérée comme une estimation de utt 


Nom 

Coefficient 

Ecart-type 

t.exp. 

Signif. % 

b 0 

160.6 

3.1 

51.40 

a 

b i 

■19,4 

3.1 

-6.20 

10,6 

b 2 

-70,6 

3.1 

-22.60 

■k 

b 3 

6,9 

3.1 

2.20 

27.9 

b 12 

■0.6 

3.1 

-0.20 

86.3 

b 13 

■3,1 

3.1 

-1.00 

50.0 

b 23 

-6,9 

3.1 

-2.20 

27.9 
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Durée 




24 h 

m 
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Conclusion 

Grandes étapes d'une étude 

• Recherche des variables importantes 

> Parmi toutes les variables susceptibles d 'influencer le phénomène étudié : 

- lesquelles ont une influence significative ou non? 

- que vaut cette influence? (positif ou négatif ) 

- y a-t-il des interactions entre les variables expérimentales ? 

• Modélisation 

> Chercher l'allure de cette influence : linéaire ou courbe; quelle équation 
permettant de décrire les variations du phénomène étudié avec les variables 
importantes 

• Optimisation 

> chercher les conditions expérimentales qui donnent le meilleur résultat pour un 
contrainte. 
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Chapitre VII 


Analyse en composante 
principale ACP 
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Variables explicatives 
et à expliquer séparées : 





l ue seule variable 
à expliquer à l:t fois 


©' 


Variables À expliquer 
xJW quantitatives 


uriuble explicatives 
quantitatives 


i'U'/i , I1C7 


Pr R. 


?L£ 



Méthodes descriptives 

* 














Problème de la Compression des Données 



Comment 

prendre en compte : 


Un grand nombre de variables Un nombre limité d’individus 


l’information pertinente étant noyée dans une 


masse d’informations partiellement redondantes 


et polluées par toutes sortes d’erreurs ? 



Peut-on réduire la masse des données 
sans perdre trop d’information ? 
et améliorer la lisibilité du contenu des données ? 


Pr R. BCHITOU 


153 







L’Analyse en omposantes Principales 

(ou rincipal Component Analysis) 


Elle consiste à transformer les Itl variables quantitatives 
initiales toutes plus ou moins corrélées entre elles 

en m nouvelles variables quantitatives non 
corrélées appelées composantes principales. 

Puis, si possible, ne conserver qu’un nombre 
k réduit de composantes. 
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Principe de l’ACP 


A partir de la matrice X des données de départ on 
construit la matrice M de corrélation : 

M = 1/n (X X) 

Avec X’ matrice transposée de X ; 

N. B. : X peut être préalablement transformée en variables 
centrées , en variables centrées réduites . 


La matrice M de corrélation est diagonalisée pour fournir 


une matrice diagonale D des valeurs propres . . À, m 


une matrice P des vecteurs propres. 
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Les vecteurs colonnes de la matrice P 
représentent les combinaisons linéaires 
des variables de départ conduisant aux 
nouveaux axes factoriels : 


les Composantes Principales. 



Les valeurs propres À, 2 v*^ m 
représentent les variances des 
individus sur les composantes 
principales correspondantes. 


D = P -1 M P 


1 


k 

1 


k 


1 


k 


Matrice P des 
vecteurs propres 


1 


k 


x 2 

Matrice^ D des 
valeurs propres 
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Le rapport de chaque valeur propre à la somme de 
toutes les valeurs propres donne la part de toute 
l'information initiale visible sur chaque axe : 

% information = 



On peut alors faire une "projection" d'un espace à k 
dimensions initiales à un espace plus réduit (2, 3) en 
sachant quel pourcentage de l'information est conservé 
(90 % par exemple). 
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Li 

L, 

- 2.124 

0.101 

- 1.420 

- 0.274 

- 0.569 

0.249 

- 0.024 

- 0.057 

0.570 

- 0.084 

1.417 

0.174 

2.149 

- 0.109 



l 2 

O 

_ c 

H 

O 

O 

1 

i 

“O 

O 

1 1 lo 1 1 ^1 


Les descripteurs dans ce nouveau système 
d’axes sont totalement décorrélés : 


Matrice de 
corrélation 



Li 

L-, 

Li 

1 


l 2 

0.00 

1 
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O 

O 


O 


O 


O 

O 


O 



Matrice de corrélation 


Soit 7 individus décrits par 



X 1 

X 2 

deux variables X. et 

x i 



1 L 

fortement corrélées 


1 



X, 

0.97 

1 
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Composantes Principales dans le cas de 2 variables 






Les données sont centrées 
autour de la moyenne 

K 

Ci = (Xi - x) 
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CR ! 

CR, 

-1.32 

-1.46 

-1.11 

-0.75 

-0.21 

-0.54 

-0.05 

0.02 

0.32 

0.43 

1.04 

0.81 

1.34 

1.48 


CR,| 

m» 

O 


O 


O 

c 


O ° 

CRj 

O 



Ou bien, les données sont 
centrées réduites 

crj = (xj - x)/a 
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Retenons le passage aux variables centrées réduites 



La seconde composante 
est orthogonale à la 
première. Pr 






Retenons le passage aux variables centrées réduites 
qui conduit à des nombres sans dimensions : 

Ce changement de variable confère le 
même poids à chacun des descripteurs. 



Le passage aux variables 
centrées réduites étant une 
transformation linéaire, il 
n’v a pas de changement de 
la matrice de corrélation. 


Matrice de corrélation 
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Valeur propre 

1.971 

0.029 

% variance 

98.571 

1.429 

% cumulé 

98.571 

100 


Vecteurs propres 


vi 

ail 

y 2 

al2 

a21 

a22 


VI 

0.707 

V2 

0,707 

0.707 

-0,707 


Calcul des coordonnées en composantes principales 

Ll = CRl*all+CR2*a21 
L2= CRl*al2 + CR2*a22 
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Les coordonnées des 
individus dans l'espace 
des composantes sont 

appelées les « SCOres ». 


L x 

l 2 

- 2.124 

0.101 

- 1.420 

- 0.274 

- 0.569 

0.249 

- 0.024 

- 0.057 

0.570 

- 0.084 

1.417 

0.174 

2.149 

- 0.109 
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Individus (axes Fl et F2 : 100.00 %) 


2.5 
2 

1.5 

? 1 

o N 

2 0.5 

— 0 

CN 

o -0.5 
x 

? - 1 
-1.5 

-2 
-2.5 

-2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 

-axe Fl (98.57 %)-> 


• i *3 ' 

.. - c *6 

i ■ i i i ■ 

• 2 

> 4 • 5 1 ^trr 


Fi 

f 2 

- 2.124 

0.101 

- 1.420 

- 0.274 

- 0.569 

0.249 

- 0.024 

- 0.057 

0.570 

- 0.084 

1.417 

0.174 

2.149 

- 0.109 
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Li 

l 2 

- 2.124 

0.101 

- 1.420 

- 0.274 

- 0.569 

0.249 

- 0.024 

- 0.057 

0.570 

- 0.084 

1.417 

0.174 

2.149 

- 0.109 



l 2 

O 

_ G 

O O . T 

i 

"O 

û 

. 

O 

1 


Les descripteurs dans ce nouveau système 
d’axes sont totalement décorrélés : 


Matrice de 
corrélation 



Li 

L, 

Li 

1 


l 2 

0.00 

1 
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A elles seules, les projections des individus sur la 
première composante : 


F — i — i *| • — i*— i — i* i *1 — - 1 i 



« visualisent ou expliquent » la quasi totalité de 
la « dispersion » des individus. 
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c: 


Le fait de négliger la part de « dispersion » expliquée 
par les projections des individus sur la seconde 
composante : 


U 


constitue une perte d’information négligeable. 
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L'information apportée par une projection, 
dépend fortement de la direction, donc ici de la 
composante principale considérée ... 
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Application 


Le levain remplace la levure boulangère utilisée dans les 
pains industriels et traditionnels. Les levures et les 
bactéries lactiques du levain font monter la pâte et 
donnent du goût au pain. 

il s'agit d'un mélange de farine , d'eau et de sel qu’on laisse 
pétrir pendant une vingtaine de minutes dans un pétrin 
mécanique pour obtenir une pâte élastique qui sera 
abandonnée une quinzaine de minutes dans le pétrin. 
Puis on a mesuré l’acidité , la viscosité, la durée et la 
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température . 





Acidité 1 

Viscosité 

Dur 

ée Tempéï 

•ature 

1 

121,5 

8100 

1 

8 

L k 

2 

3 

126 

124 

9360 

13000 

1 

1 

8 

22 

L2X 

L1 \ 

4 

131,5 

12260 

1 

22 

L2 

5 

147,5 

6260 

7 

8 

L1 

6 

149,5 

8750 

7 

8 

L2 

7 

153,5 

11250 

7 

22 

L1 

8 

154,25 

12766 

7 

22 

L2 

9 

149,5 

7633 

15 

8 

L1 

10 

150,5 

9366 

15 

8 

L2 

11 

155 

8800 

15 

22 

L1 

12 

157,75 

11100 

15 

22 

L2 

13 

152 

6800 

21 

8 

L1 

14 

155 

9466 

21 

8 

L2 

15 

157 

94<% R 

. bAhitou “ 

U 

16 

159 

9700 

21 

22 

L2 
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CORRELATIONS DES VARIABLES INITIALES 

ACIDITE VISCOSITE CONSERVATION TEMPERATURE 


ACIDITE 

1 




VISCOSITE 

-0,02 

1 



CONSERVATION 

-0,072 

-0,341 

1 


TEMPERATURE 

-0,194 

0,726 

0 

1 


ANALYSE DES VARIABLES SYNTHETIQUES 



AXE 1 

AXE 2 

AXE 3 

AXE 4 

ACIDITE 

0,2 

0,73 

-0,65 

0,05 

VISCOSITE 

-0,94 

0,14 

-0,12 

-0,3 

CONSERVATION 

0,38 

■0,71 

-0,58 

-0,12 

TEMPERATURE 

-0,87 

■0,29 -0,27 
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AXE 2 
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Indiv. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 


Axe 1 
0,74 
0,32 
- 2,2 
-1,91 
1,69 
0,84 
-1,27 
-1,79 
1,51 
0,91 
-0,13 
-0,91 
2,02 
1,11 
-0,71 
- 0,21 


Axe 2 

0,73 

0,89 

0,57 

0,66 

0,57 

0,77 

0,48 

0,59 

0,03 

0,16 

-0,33 

-0,13 

-0,49 

-0,26 

-3,55 

-0,7 


Axe 3 

1,47 

1,31 

0,54 

0,45 

0,62 

0,43 

-0,38 

-0,49 

-0,14 

-0,27 

-0,89 

-1,09 

-0,62 

-0,84 

1,42 

-1,51 


Axe 4 
0,22 
- 0,2 
- 0,22 
-0,06 
0,72 
-0,14 
0,26 
-0,27 
-0,04 
-0,64 
0,82 
0,03 
0,04 
- 0,88 
-0,06 
0,3 
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Axe 2 


® 4 * 8 n 

2 * ®* 6 ® 
1 # 

12 © 
♦ 

1 1 

16 




Axe 1 
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